Near-field heat transfer in a scanning thermal microscope 
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We present measurements of the near-field heat transfer between the tip of a thermal profiler 
and planar material surfaces under ultrahigh vacuum conditions. For tip-sample distances below 
10 -8 m our results differ markedly from the prediction of fluctuating electrodynamics. We argue 
that these differences are due to the existence of a material-dependent small length scale below 
which the macroscopic description of the dielectric properties fails, and discuss a corresponding 
model which yields fair agreement with the available data. These results are of importance for the 
quantitative interpretation of signals obtained by scanning thermal microscopes capable of detecting 
local temperature variations on surfaces. 
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Radiative heat transfer between macroscopic bodies in- 
creases strongly when their spacing is made smaller than 
the dominant wavelength Ath of thermal radiation. This 
effect, caused by evanescent electromagnetic fields exist- 
ing close to the surface of the bodies, has been studied 
theoretically already in 1971 by Polder and van Hove 
for the model of two infinitely extended, planar sur- 
faces separated by a vacuum gap [l[, and re- investigated 
later by Loomis and Maris [2| and Volokitin and Pers- 
son 0, |5J ■ While early pioneering measurements with fiat 
chromium bodies had to remain restricted to gap widths 
above 1 pi |H, and later studies employing an indium 
needle in close proximity to a planar thermocouple re- 
mained inconclusive Q, an unambiguous demonstration 
of near-field heat transfer under ultrahigh vacuum con- 
ditions and, thus, in the absence of disturbing moisture 
films covering the surfaces, could be given in Ref. [7|. 

The theoretical treatment of radiative near-field heat 
transfer is based on fluctuating electrodynamics §■ 
Within this framework, the macroscopic Maxwell equa- 
tions are augmented by fluctuating currents inside each 
body, constituting stochastic sources of the electric and 
magnetic fields E and H. The individual frequency com- 
ponents j(r,w) of these currents are considered as Gaus- 
sian stochastic variables. According to the fluctuation- 
dissipation theorem, their correlation function reads 

(j Q (r,^(r',a/)) 

= -E(u, p) e" (oj) 5 a0 5(r - r') 5(u> - u/) , (1) 

where E(u, (3) — hid / (exp((3hu>) — l), with the usual in- 
verse temperature variable /3 — l/(fceT); the angular 
brackets indicate an ensemble average. Moreover, e"(u>) 
denotes the imaginary part of the complex dielectric func- 
tion e(cj) = e'(ui) + ie"(cj). It describes the dissipative 
properties of the material under consideration, which is 
assumed to be homogeneous and non-magnetic. Thus, 
Eq. ([TJ) contains the idealization that stochastic sources 
residing at different points r, r' are uncorrelated, no mat- 
ter how small their distance may be. Applied to a mate- 



rial occupying the half-space z < 0, facing the vacuum in 
the complementary half-space z > 0, these propositions 
can be evaluated to yield the electromagnetic energy den- 
sity in the distance z above the surface, giving [lOj 
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Here, the densities pb and pn symbolically specify the 
electric and magnetic contribution, respectively; ru de- 
notes the Fresnel amplitude reflection coefficient for TM- 
modes with wave vector of magnitude ujk/c parallel to the 
surface. The symbol X± abbreviates the corresponding 
terms for TE-modes. The wave vector oriented normal to 
the surface, of magnitude Ljp/c, distinguishes propagat- 
ing modes with real p = \/l — k 2 for k < 1 from evanes- 
cent modes with imaginary p = is/ k 2 — 1 for k > 1. 

Expression ([2]) for the energy density, obtained strictly 
within the framework of macroscopic electrodynamics, 
diverges for small distances z from the surface; for 
z/Xth <C 1) onc finds the power law (u(z)) cx z~ 3 
Hence, it has been suggested that also the energy dis- 
sipated in the tip of a tiny probe close to the surface 
should scale inve rsely proportional to the cube of the tip- 
sample distance [111 |l2j. However, the entailing diver- 
gence clearly is not borne out by the actual physics fl3l — 
Il5j | . The divergence may formally be avoided by replacing 
the upper boundary of integration, k = oo in Eq. (J2J) , by a 
finite cutoff k c , thereby excluding the problematic large- 
wave number contributions to the "evanescent" part of 
the energy density It is important to note that the 
divergence of the energy density (O close to the material 
surface reflects a shortcoming of the underlying macro- 
scopic theory: Considering a metal, the dielectric prop- 
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erties of which are largely determined by the conduction 
electrons, one expects that any contributions from spatial 
Fourier components shorter than their mean free path are 
inadequately dealt with [lj]. More generally, the spatial 
delta-like correlation (TjQ) becomes problematic on length 
scales such that the microscopic properties of the mate- 
rials start to make themselves felt. These observations, 
in their turn, imply that experiments on fluctuating elec- 
tromagnetic fields in the extreme near-field regime, where 
traditional macroscopic fluctuating electrodynamics can 
no longer be taken for granted, may yield important in- 
formation on microscopic material properties. 

In this Letter, we report on measurements of the near- 
field heat transfer between the tip of a scanning thermal 
microscope and surfaces of gold (Au) or gallium nitride 
(GaN). We have fabricated a thermosensor, integrated 
into the tip of a variable-temperature scanning tunneling 
microscope (VT-STM) , which allows us to determine the 
heat transfer even for tip-sample distances on the order 
of I nm. We argue that our sensor essentially probes the 
near-field energy density close to the sample, and demon- 
strate that the experimental data differ markedly from 
the standard prediction ([2]) for distances below 10 nm. 
A simple, but physically motivated ansatz for the de- 
scription of the short-range dielectric material properties 
then leads to qualitative agreement with the measured 
data, allowing one to extract material-dependent length 
scales L below which the macroscopic theory fails. 

When assessing the near-field heat flux between two 
bodies of different temperature, precise location of their 
positions of zero separation is of key importance. Since 
no body has a mathematically flat surface, this is to some 
extent a matter of definition. We have chosen to record 
the heat transfer between a cooled sample and the tip 
of a VT-STM at nearly room temperature, so that zero 
separation of the two surfaces corresponds to a certain 
level of electron coupling, i.e., to a certain tunnel cur- 
rent. To exclude any mechanism of heat transfer other 
than radiation, one has to work under ultrahigh vacuum 
(UHV) conditions. Otherwise, any surface adsorbate, or 
surrounding gas, might result in additional contributions 
to the heat transfer, masking the radiative effect. 

The heat flux between the warm tip and the cooled 
sample is measured through the resulting slight diminu- 
ation of the temperature of the very tip compared to the 
rest of the sensor. Since small temperature differences 
have to be detected over a small sensor, any self- heating 
has to be carefully avoided. Therefore, we employ a ther- 
mocouple integrated into the tip of our VT-STM. As 
sketched in Fig. [T] (a), a thin platinum wire has been 
melted into a glass micropipette. Subsequently, the part 
of the wire protruding from the pipette has been electro- 
chemically etched to form a sharp tip [lf|. The pipette 
has then been covered by a gold film with a thickness of 
about 25 nm, having electrical contact with the platinum 
wire only at the very end of the tip. This end thus forms 
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FIG. 1: (a) Cross section of the micropipette glued into a 
tip holder. The thermoelectric voltage Vth builds up between 
the inner platinum wire and the outer gold film. The tunnel 
potential is applied between the sample and the grounded gold 
film, (b) Dependence of the thermovoltage on the absorbed 
power AP of calibrating laser light for two different sensors. 



the measuring contact of the resulting coaxial thermo- 
couple, while the reference contact is located in the back 
at the support of the micropipette, with good thermal 
coupling to the surrounding which acts as a heat bath. 

The support of the micropipette is sitting in the scan- 
ner of a commercial VT-STM in an UHV chamber made 
by Omicron, held by a ring-shaped magnet. In the hole 
of the magnet a gold-plated, spring-loaded contact pin 
makes electrical contact with the platinum wire, thus 
forming the reference contact. The gold film is grounded, 
whereas the sample is at tunnel potential, in order to de- 
couple the tunneling signal from the thermoelectric volt- 
age Vth- This voltage is first amplified by a Keithley 
nanovolt preamplifier (model 1801) and then measured 
by a high resolution multimeter (Keithley model 2001). 
The temperature of the sample is lowered during the mea- 
surements with liquid nitrogen via a coldfinger to 100 K, 
establishing a temperature difference between the tip and 
the sample surface of about 200 K. 

The Seebeck coefficient S = Vth/ AT of our sensor, 
quantifying the ratio of the generated thermoelectric 
voltage and the temperature difference AT between the 
two contacts of the thermocouple, is determined with a 
setup consisting of a droplet of oil held by a small heat- 
ing coil made of tungsten wire. The temperature of the 
oil is measured by a commercial type-K thermocouple 
reaching from one side into the droplet, while the sensor 
enters it from the other side. The oil temperature can 
be varied by changing the current through the coil. We 
obtain S = 8 uV/K at room temperature, which is close 
to the value found in literature [171 ] . 

The heat resistance R& = AT/ AP of the sensor, re- 
lating the heat power AP absorbed or emitted by the 
tip to the resulting temperature difference between the 
contacts, was determined by placing the tip in the focus 
of a 1 mW cw laser diode (wavelength 670 nm). The 
fraction of the light power which did not hit the tip 
surface was measured by a power meter positioned be- 
hind it. The absorbed power was estimated, according 
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FIG. 2: Measured heat current AP (in Watts) between the 
microscope tip and a gold layer (circles) vs. tip-sample dis- 
tance z. The dashed line, which coincides with the full one 
for larger z, corresponds to the prediction AP t h of standard 
fluctuating electrodynamics, based on Eq. p|). The full line 
is obtained from Eq. <(Sj with the modified dielectric func- 
tion ©, setting L = L t i P = 1.2 ■ KT 8 m. 



to AP = (P - P) • (1 - R), from the difference P - P 
of the power recorded without and with the tip being 
present, and the reflectivity R = 0.96 of its gold surface 
at 670 nm 0. The expected linear dependence of the 
thermovoltage on the absorbed power is well confirmed 
in Fig. [T] (b), showing our results for two different tips. 
From the slopes 0.18 pV/pW and 0.43 pV/pW we obtain 
heat resistances of 23 K/ mW and 54 K /mW, respectively. 
Knowing both a sensor's Seebeck coefficient S and its 
heat resistance Pth, one can deduce the near-field heat 
flux AP between the tip and a closely spaced sample of 
different temperature from the observed thermovoltage, 
according to AP = V th /(S'P th ). Measurements of the 
distance dependence of the heat transfer were performed 
by retracting the STM tip from the tunnel distance, while 
the distance itself was determined by means of the piezo 
coefficient of the scanner. Results of such measurements 
are depicted in Fig. [2] for a sample consisting of a gold 
layer, and in Fig.|3]for a sample of GaN. In both cases, the 
sensor with Pth = 54 K/mW has been employed. Dur- 
ing these measurements, we have carefully checked that 
the crosstalk between the tunnel current signal and the 
thermovoltage remains negligibly small. The absence of 
interference is indicated by the fact that the tunnel cur- 
rent decreases strongly in a range of distances where the 
observed thermovoltage stays almost constant. 

A theoretical discussion of the heat transfer between 
an idealized tip and a flat surface, which may serve as 
a guideline for the analysis of our data, has been given 
by Mulet et al. [l2| ■ These authors have modeled the tip 
by a small dielectric sphere of radius r and assumed the 
incident electric field to be uniform inside the sphere, so 
that it acts like a point-like dipole. If the temperature of 
the sample is significantly lower than that of the tip, as in 
our case, the heat current flowing back from the sample 
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FIG. 3: As Fig. [2] for a sample of GaN, setting L = 1.0 • 
10~ 10 m and L tip = 1.2 • 10~ 8 m. 

to the tip can be neglected. The total flux between the 
surface and the tip then is determined entirely by the 
current directed from the tip to the sample, according to 

/>oo />OC 

AP = duj dn cv(w)pe(w, k, (3, z) , (3) 
Jo Jo 

where a(uj) — 2w(47rr 3 )e" ip /|e t ip+2| 2 describes the dielec- 
tric properties of the sphere, and the temperature enter- 
ing ,0e is that of the tip. Taking this expression at a rep- 
resentative frequency uj , one has AP « a(o;o)eo(E 2 )/2 , 
so that, within the scope of the model, the heat flux reg- 
istered by the tip should be proportional to the electrical 
energy density of the flat sample, evaluated, however, at 
the temperature of the tip. 

For distances larger than about 10~ 8 m, our exper- 
imentally observed heat transfer is, to good accuracy, 
proportional to the total energy density as given by 
Eq. ([2]), not to the electric field contribution alone. 
Since the constant of proportionality, which carries the 
dimension of area times velocity, may differ substan- 
tially from a(wo), we focus on the scaled energy density 
AP t h := ira 2 c(u(z)}, where c is the velocity of light, and 
employ the effective sensor area ira 2 as a fitting param- 
eter. Modeling the dielectric function e(ui) for Au by a 
Drude ansatz with parameters taken from Ref. [l||, and 
that for GaN by the "reststrahlen" -formula with param- 
eters from Ref. j20j, we obtain the dashed lines in Figs. [5] 
and [3l setting a = 60 nm. This value is in accordance 
with scanning electron microscopy studies of the tip, and 
describes both experimental data sets for z > 10 nm, as 
it should. The latter fact also indicates that the use of 
Eq. ©, i.e., the neglect of the field's distortion by the 
tip, is justified here. 

In the case of GaN, the theoretical curve for AP t h di- 
verges as z~ 3 for sensor-sample distances below 10 nm. 
In contrast, for Au this familiar behaviour would become 
apparent only at substantially smaller z Q. However, the 
experimental data clearly show a different trend, leveling 
off to values which for the smallest accessible distances 
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are significantly lower than AP t h- We interpret this find- 
ing as evidence for the short-distance deficiency of the 
macroscopic theory, as expressed by the delta-like corre- 
lation function (JXJ) of the stochastic source currents: In a 
real sample, there is some finite correlation length L. 

In principle, one should then also account for non-local 
effects, which requires distinguishing a transversal and a 
longitudinal part of the permittivity 2l(. Instead, here 
we propose a more simple, qualitative approach: On the 
one hand, only the imaginary part e"(w) of the response 
function e(cu) enters into the fluctuation-dissipation the- 
orem and thus into the correlation function ([IJ; on the 
other, the Kramers-Kronig formula relates the imaginary 
part to the real one. Hence, a plausible and consistent 
ansatz for an effective permittivity depending explicitly 
on the transversal wave number is 



In summary, we have obtained experimental data for 
the near-field heat transfer between a thermal profiler 
and flat material surfaces under UHV conditions. We 
have reached the extreme near-field regime, where the 
variation of the heat transfer rate with the distance be- 
tween microscope tip and sample differs distinctly from 
the divergent behaviour predicted by standard macro- 
scopic fluctuating electrodynamics, and have interpreted 
our observations in terms of finite microscopic correla- 
tions inside the materials. While the shortcomings of the 
macroscopic theory are, in principle, well known [1,0], 
their manifestation in an actual experiment indicates a 
still unexplored potential of thermal microscopy as a new, 
quantitative tool for the nanometer-scale investigation of 
solids. 



e(w, k) := 1 + [e'M - 1] /(«) + ie»/(/c) , (4) 

where the function J(k) accounts for the lateral corre- 
lations, such that it approaches unity and thereby re- 
stores the local case when ujk/c <C L~ x , but vanishes for 
large wave numbers, when ujk/c 3> L . As a convenient 
guess, we take a Gaussian f(n) = exp ( — (Lluk/c) 2 ) , and 
consider L as a parameter to be determined by fitting the 
data. This parametrization (j4]) has the distinct advan- 
tage that the Maxwell equations for systems with plane 
translational invariance [jj remain formally unchanged; 
it is only that e(w) has to be replaced by e"(u;, k). In par- 
ticular, the energy density can again be obtained from 
Eq. ([2]), if only the reflection coefficients rn and r± are 
adapted in this manner. 

Besides the dielectric properties of the sample, also 
those of the sensor enter into the data, as exemplified 
by the dipole model ([3]). Hence, we have to introduce 
both a correlation length L of the sample and a further 
correlation length L t ip of the sensor, and parametrize the 
experimentally observed heat current in the form 



AP 



(5) 

Using this ansatz ([5]), we finally obtain the full curves in 
Figs. H and El setting L = 1.2 • 1CT 8 m for Au and L = 
1.0 • 1CT 10 m for GaN, while L tip = 1.2 • 10" 8 m in both 
cases. These curves capture the experimental data quite 
well, thus lending strong support to our line of reasoning. 
It is also encouraging to observe that the numerical value 
of L obtained for Au indeed turns out to be on the order 
of the mean free path of electrons in metals, whereas that 
for GaN is considerably shorter, as it should. Although 
the thermally relevant component of our sensor probably 
is confined to the Au layer, its correlation length not 
necessarily has to coincide with that of the gold sample, 
as it actually does in our case, but might be geometrically 
restricted in alternative setups. 
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